Energy conservation for point particles undergoing radiation 

reaction 



Theodore C. Quinn and Robert M. Wald 
Enrico Fermi Institute and Department of Physics 
University of Chicago 
5540 S. Ellis Avenue 
Chicago, Illinois 60637-1433 
(February 7, 2008) 



Abstract 



For smooth solutions to Maxweh's equations sourced by a smooth charge- 
current distribution ja in stationary, asymptoticahy flat spacetimes, one can 
^ I prove an energy conservation theorem which asserts the vanishing of the sum 

of (i) the difference between the final and initial electromagnetic self-energy 
' of the charge distribution, (ii) the net electromagnetic energy radiated to in- 

■ finity (and/or into a black hole/white hole), and (iii) the total work done by 

^\ . the electromagnetic field on the charge distribution via the Lorentz force. A 

similar conservation theorem can be proven for linearized gravitational fields 
off of a stationary, asymptotically fiat background, with the second order 
Einstein tensor playing the role of an effective stress-energy tensor of the lin- 
W)' earized field. In this paper, we prove the above theorems for smooth sources 

and then investigate the extent to which they continue to hold for point par- 
^ ] tide sources. The "self-energy" of point particles is ill defined, but in the 

H ' electromagnetic case, we can consider situations where, initially and finally, 

the point charges are stationary and in the same spatial position, so that 
the self-energy terms should cancel. Under certain assumptions concerning 
the decay behavior of source-free solutions to Maxwell's equations, we prove 
the vanishing of the sum of the net energy radiated to infinity and the net 
work done on the particle by the DeWitt-Brehme radiation reaction force. 
As a by-product of this analysis, we provide a definition of the "renormalized 
self-energy" of a stationary point charge in a stationary spacetime. We also 
obtain a similar conservation theorem for angular momentum in an axisym- 
metric spacetime. In the gravitational case, we argue that similar conservation 
results should hold for freely falling point masses whose orbits begin and end 
at infinity. This provides justification for the use of energy and angular mo- 
mentum conservation to compute the decay of orbits due to radiation reaction. 
For completeness, the corresponding conservation theorems for the case of a 
scalar field are given in an appendix. 
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I. INTRODUCTION 



The problem of calculating the motion of an isolated body coupled to fields in curved 
spacetime is an old one which is currently receiving renewed interest. This recent interest 
is largely driven by the need for accurate calculations of processes which emit gravitational 
waves (in anticipation of results from the new generation of gravitational wave detectors) 
but similar issues arise for bodies coupled to electromagnetic and scalar fields. In order to 
describe in a simple manner those aspects of a body's motion which are independent of its 
detailed internal structure, one often attempts to calculate the motion of a "point particle." 
A central problem in all such investigations is to calculate the effects of the particle's own 
fields, commonly referred to as "self force" or "radiation reaction" effects. This problem is 
mathematically ill posed since the fields diverge on the world line of the particle itself. 

Nevertheless, there is a long history of attempts to calculate the self-force on a particle 
directly from the local fields. All such schemes involve some prescription for subtracting 
away the infinite contributions to the force due to the singular nature of the field on the 
particle's world line. In 1938, Dirac produced a force expression for a point charge coupled 
to an electromagnetic field in Minkowski spacetime by imposing local energy conservation 
on a tube surrounding the particle's world line The infinite contributions to the force 
were subtracted through a "mass renormalization" scheme. In 1960, Dewitt and Brehme 
generalized this approach to an arbitrary curved background spacetime. (A trivial calcula- 
tional error in their paper was later corrected by Hobbs 0.) More recently, Mino et al. |^ 
further adapted this approach to produce an expression for the self force on a massive parti- 
cle coupled to linearized gravity on a vacuum background spacetime. Recently, we also have 
derived |^ the formulas for the electromagnetic and gravitational self-forces by using an 
axiomatic approach which, in effect, regularizes the forces by comparing forces in different 
spacetimes. 

Despite the fact that it is thereby known, in principle, how to calculate the electromag- 
netic and gravitational self-force on a point particle, serious difficulties arise in practice when 
one attempts to evaluate this self-force on account of the difficulties in computing the "tail 
term" contribution. Indeed, the evaluation of the "tail term" is highly nontrivial even in 
the slow motion, weak field limit PJ^. Consequently, many researchers have employed the 
following iterative strategy for calculating the motion of point particles. First, one calcu- 
lates the motion of the particle in the absence of self- force effects. Then, the energy (and/or 
angular momentum) radiated to infinity by the resulting fields is calculated. Finally, the ef- 
fects of this energy (and/or angular momentum) loss are introduced as a perturbation to the 
particle's motion. Note that the applicability of this approach is limited in that the energy 
and angular momentum of the particle are not even defined when the background spacetime 
fails, respectively, to be stationary and axisymmetric. In addition, even when the energy 
and/or angular momentum of the particle are well defined, the energy and angular momen- 
tum radiated to infinity would be expected to equal the energy and angular momentum loss 
by the particle only in a time averaged senseQ Hence, the approach is, in essence, limited to 



The example of a point charge in Minkowski spacetime which undergoes a period of uniform 
acceleration explicitly demonstrates the failure of temporally local conservation of this sort, since 
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calculating the secular decay of an otherwise stable orbit. Furthermore, in the absence of 
spherical symmetry, the loss of energy and angular momentum in general is not sufficient 
to determine even this secular decay (for example, it does not determine the variation of 
the Carter constant for non-equatorial orbits in the Kerr spacetime). Finally, this approach 
does not allow one to calculate any so-called "conservative forces", i.e., contributions to the 
self-force which are not associated with energy or angular momentum loss. Nevertheless, 
this approach is extremely simple to apply and has been widely used to estimate the effects 
of radiation reaction on the motion of point particles. 

Although the derivations of the self- force given in |]T|, and ||^ were heuristically 
motivated by local conservation of energy, it is far from obvious, a priori, that they satisfy 
the property of "global energy conservation" as assumed in the iterative procedure described 
above. Specifically, in order to justify the iterative procedure, it is necessary that the total 
energy and angular momentum radiated to infinity by the electromagnetic or gravitational 
fields be equal, respectively, to the net work and torque done by the self-force over the world 
line of the particle. The main purpose of this paper is to investigate the extent to which 
this is the case. 

Our main results are the following: In the electromagnetic case, we consider stationary, 
globally hyperbolic, asymptotically fiat spacetimes. We assume that source-free Maxwell 
fields satisfy a certain decay property, and that the advanced and retarded solutions with 
stationary sources have suitable fall-off at infinity. We then prove that if a point charge 
is asymptotically stationary in the past and future and in the same position, then the net 
electromagnetic energy radiated to infinity (and/or into a black hole/white hole) is equal 
to the integral over the particle's world line of the force expression given by DeWitt and 
Brehme (as corrected by Hobbs) contracted with the timelike Killing field. This provides 
some justification for the use of the above global energy conservation method to calculate 
the effects of radiation reaction. (Alternatively, the fact that the DeWitt-Brehme force gives 
rise to global energy conservation may be viewed as providing further justification for its 
own validity.) As byproduct of this theorem, we show that it is possible to consistently 
define the "renormalized self-energy" of a stationary point charge. We argue that a similar 
energy conservation result should also hold for a point charge which is in inertial motion 
near infinity in the asymptotic past and future. We also show that similar results hold for 
angular momentum conservation in axisymmetric spacetimes. As an additional byproduct of 
our analysis, we show that our conservation theorems also hold for the force expression used 
by Gal'tsov and others, in which the Lorentz force associated with the advanced-minus- 
retarded "radiative" field is used. This result clarifies the relationship between the DeWitt- 
Brehme and Gal'tsov formulas and demonstrates explicitly that global energy conservation 
alone is insufficient to determine the local force. 

Our analysis of the gravitational case closely parallels that of the electromagnetic case. 



the radiation reaction force (and, hence, the work done by it) vanishes at retarded times during 
which a nonzero flux of radiation reaches infinity. Nevertheless, the net work done by the radiation 
reaction force equals the net energy radiated to infinity if the motion of the particle is static at early 
and late times [^. The generalizations of such average conservation results to curved spacetime 
and to the gravitational case are the main subjects of the present paper. 
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with the second-order Einstein tensor |T0| of the metric perturbation replacing the stress- 
energy tensor of the electromagnetic field. However, our analysis is hampered by a number 
of technical obstacles — such as the fact that the energy flux as calculated from the second 
order Einstein tensor has been proven to agree with the Bondi flux only for perturbations 
of compact spatial support — and its domain of applicability is limited to freely falling point 
masses (since if the point mass is not freely falling, an additional stress-energy source for 
the metric perturbation must be present). Nevertheless, we argue that energy and angular 
momentum conservation results analogous to the those in the electromagnetic case should 
hold for a point mass which is in geodesic motion near infinity in the asymptotic past and 
future. 

The electromagnetic case is analyzed in Sec. II. We begin by proving energy conservation 
theorems for the case of smooth charge-current distributions in a spacetime with no black 
or white holes. The desired theorem for the case of point particle sources is then proven. 
Generalizations to establish angular momentum conservation and to allow for the presence 
of black and white holes are then described. 

Linearized gravitational perturbations are considered in Sec. III. Again, we begin by 
proving a conservation theorem for smooth stress-energy sources and then analyze the point 
particle case. 

Finally, for completeness, we present the analogous results for a scalar field in an ap- 
pendix. 



Our notation and conventions throughout the paper follow [11 



II. ELECTROMAGNETIC CASE 

In this section, we will state and prove our conservation theorems for the case of a Maxwell 
field in a fixed background spacetime. For simplicity, we shall first consider conservation 
of energy in a stationary, asymptotically fiat spacetime which contains no black or white 
holes. Generalizations to other symmetries and to allow for the presence of black holes and 
white holes will be discussed at the end of this section. Thus, until these generalizations are 
considered, we will restrict consideration to spacetimes satisfying the following conditions: 
Spacetime assumptions: Let (M , gab) be a globally hyperbolic spacetime that is asymp- 



totically flat at null and spatial infinity in the sense of Ashtekar and Hansen / [7^ and is 
stationary in the sense that it possesses a Killing field that is asymptotically a time trans- 
lation at infinity. We further assume that no black holes or white holes are present in 
{M,gab) (i.e., the domain of outer communications is the entire spacetime) and that there 
exists a smooth, spacelike Cauchy surface, S for (M, gab) such that in the unphysical space- 
time {M,gab), S U 2*^ is compact. (This implies that S is of the form of a disjoint union 
S = Send U S' where Send has the topology ofTZ^ minus a closed ball and S' is compact.) 

Without loss of generality, we may assume that f^Ua < everywhere on E, where n° is the 
future directed normal to S (see proposition 4.1 of [0). We now deform E to the future in a 
small neighborhood of i° so that the deformed surface, C"^, satisfies [C^HE] D E' and, in the 
unphysical spacetime, C+ remains a smooth, spacelike hypersurface, but now intersects X+ 
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in a cross-section ^ We similarly deform S to the past to construct a smooth, spacelike 
hypersurface C~ such that [C~ flS] D E' and C~ intersects X~ in a cross-section S~ . For all 
t"*" > 0, we define C~^{t~^) to be the hypersurface obtained by "time translating" along 
the orbits of t° by t~^. Similarly, for all t~ < 0, we define C^(t^) to be the hypersurface 
obtained by "time translating" C~ along the orbits of by . It then follows that for all 
t+>0 and t- < the surface C+(t+) U [J+ n J-{S+{t+))]UC-{t-) U [J- n J+{S-{t~))]Ui^ 
bounds a compact region V{t~^,t^), of the unphysical spacetime, as illustrated in Figure |I|. 
Furthermore, it follows that any p & M lies within V{t~^,t~) for sufficiently large t~^,t~. 

In a spacetime (M, gab) satisfying the above properties, we wish to consider solutions to 
Maxwell's equations with source ja, 

W'^Fab = -47rj, (1) 

"^[aFbc] = 0. (2) 

It follows immediately from Maxwell's equations that the stress-energy tensor 

Tab = (^Fa'^Fbc — -^dabF'^'^Fcdj (3) 

satisfies 

V'Tab = -Fabf. (4) 

We define the stress-energy current three-form Jabc by 

Jabc = t'^T^e^abc^ (5) 

where eabcd is the (positively oriented) volume element associated with the spacetime metric. 
We choose the orientation of all spacehke and null hypersurfaces to be given by v°'eabcd, where 
is any future directed timelike vector field. Then, the integral of Jabc over a spacelike 
hypersurface represents the total electromagnetic energy on that hypersurface, whereas the 
integral of Jabc over a null hypersurface represents the fiux of electromagnetic energy through 
that hypersurface. Equation @ implies that 

{dJ)abcd = —'^^{t''Tfe)eabcd 
= —f'V^Tfe'^abcd 

= fFefJ^ eabcd (6) 

Initially, we shall consider smooth solutions to Maxwell's equations with a smooth source 
ja, but later in this section we will consider solutions with point particle sources which are 
smooth away from the world line of the particle. In all cases we shall consider only solutions 
which satisfy the following three conditions: 



^One way of doing this would be to choose a cross-section 5+ of Z+ such that [J^ (5+) fl S] C Send 
and consider the C^, partially null hypersurface [S — U [J^{S^) n J+(S)]; then smooth 

this hypersurface (see [14|) to a smooth spacelike hypersurface. 
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Maxwell field assumptions: (1) ja vanishes in a neighborhood ofI+ U T" U i". (2) The 
unphysical Maxwell field Fab = Fab continuously extends to X"*" and 1~ . (3) The physical 
energy current three-form t^Tde^abc^ fo>lls off sufficiently rapidly at spatial infinity that (a) its 
integral over Send finite and (b) its integral over Br r\ M vanishes as r ^ 0, where Br is 
a suitably chosen coordinate sphere of radius r around i^ in the unphysical spacetime. 

Condition (2) lias the immediate consequence that the unphysical energy current three- 
form t'^Tde^abc'' (where labc'^ is the metric-compatible volume element associated with the 
unphysical metric, with index raised using the unphysical metric) continuously extends to 
X"*" and X~. However, for a conformally invariant field such as Fab, the energy current three- 
form is conformally invariant with conformal weight 0, so condition (2) implies that the 
physical energy current three- form f^Tde^abc^ also continuously extends to X"*" and X^. 

Our fundamental result is simply a direct application of Stokes' theorem to the integral of 
eq. (H) over the region V{t'^ ,t~). Taking account of the absence of contributions from spatial 
infinity as a consequence of assumption (3) above (i.e., initially excluding from V{t'^,t~) a 
coordinate ball of radius r around i° and then letting r — > 0), we obtain 

/ t Tabf^cde / ^ Tab(^cde 

JC+{t+) JC-{t-) 

b I ^arp , b 



+ / t"'Tab^cde — I t"'Tab€cde — / f^Fabj ^cdef-, (7) 

Ji+nj-{s+{t+)) Ji-nJ+{s+(t-)) Jv(t+,t-) 

Taking the limits as t"*" ^ oo and t^ — oo, we obtain our desired theorem for the case of 
smooth sources: 

Theorem 2.1 Let {M,gab) satisfy the properties stated above and let Fab be a smooth 
solution to Maxwell's equations with smooth source ja which satisfies the three properties 
stated above. Then, 

lim / f'Tabecde' - lim / f'Tabecde' 

t+^+oo Jc+{t+) t-^-oo Jc-{t-) 

+ / f^Tab^cde — / t^'Tab^cde = / t"" Fabj^ ^cdef i (§) 

JJ+ Ji- Jm 

provided that each of the limits and integrals in the above formula exist. 

The above theorem has a straightforward physical interpretation. The first pair of terms 
on the left side can be interpreted as being the difference between the initial and final 
electromagnetic energy of the charge-current distribution. The second pair of terms is just 
the net electromagnetic energy radiated to infinity. Finally, since Fabj^ is just the Lorentz 
force acting on the charge distribution, the right side is the negative of the net work done by 
the electromagnetic field on the charge distribution, or, equivalently, it is the net work done 
by the charge distribution on the electromagnetic field. This, in turn, has the interpretation 
of representing the net amount of "mechanical energy" which is converted to electromagnetic 
energy. Thus, eq. (H) can be interpretated as stating that total energy is conserved. 

Unfortunately, this theorem as stated does not have a direct counterpart for point particle 
sources. It is true that the net electromagnetic energy radiated to infinity is perfectly well 
defined in the limit of point particle sources. Furthermore, although the integrand on the 
right side of eq. (^ becomes singular in the point particle limit, one might nevertheless hope 
that, in a suitable limiting process, the integral on the right side of eq. (§) would converge 
to 
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Jt^fadT (9) 

where is the DeWitt and Brehme [0] expression for the electromagnetic self-force on the 
particle (see eq. (^) below). However, the self-energy terms in eq. (§) become hopelessly 
divergent in the point particle limit, so, as it stands, eq. (||) will not make sense for point 
particles. 

Nevertheless, we can proceed by considering cases where the self-energy terms in eq. 
should cancel. In particular, this should occur if at sufficiently early times the charge-current 
source ja is stationary (i.e., £tja = 0), and at sufficiently late times, ja returns to the same 
stationary state, i.e., if ja differs from a stationary distribution only in a compact spacetime 
region. However, even in this case, in order to ensure cancellation of the self-energy terms, 
it is necessary to assume suitable decay properties of source free solutions, as well as some 
properties of stationary solutions. 

It order to formulate our decay assumptions, we need to introduce a suitable norm, 
and it is most convenient to introduce this norm in terms of the unphysical variables. On 
the hypersurface C"*", we decompose the unphysical electromagnetic field Fab = Fab into its 
electric and magnetic parts with respect to the unphysical unit normal = f2~^n" to C^, 

Fa = Fabu' (10) 



Ba = -hab'^F^n' (11) 

On C^, we define 

\\F\\=^nv{E'^Ea + B'^BaY'^ (12) 

c+ 

where the indices are raised and lowered with respect to the unphysical metric. Note that for 
any Fab which is nonsingular in the physical spacetime, we have < oo provided only that 
Fab continuously extends to X"*", as we have assumed. Note also that since the components 
of t"" are bounded in an unphysical orthonormal frame associated with n"-, it follows that 
everywhere on C"^ we have \Tabt°''n}'\ < C||-F||^ for some constant C. Since has finite 
volume with respect to the unphysical metric, it follows that there exists a constant, c, such 
that 

/ f'Tablcde <C\\F\\'' (13) 

Jc+ 

Since the energy current 3-form is conformally invariant, this is equivalent to 

Jj^Tabe^de' < c\\Ff (14) 

i.e., the norm we have introduced on C'^ bounds the energy on C"^. 
In terms of the physical variables, we have 

= snpn-^iE'^Ea + B'^Baf'^ (15) 
c+ 
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where Ea and Ba are now defined with respect to the physical unit normal and indices 
are now raised and lowered with respect to the physical metric. We could similarly define a 
norm on C~^(t~^) by eqs. (^2|) or (0), but this would not be useful on account of the possible 
"time dependence" of the conformal factor Q. Thus, we instead define the function Q' on 
C^(t^) by Lie transport of Q along f^. ( Equivalent ly, we could require that the conformal 
factor Q defining the unphysical spacetime be chosen so as to satisfy £t^ = to the future 
of C+.) For all t+ > 0, we define 

||F||(t+)= snp {n'r^E'^Ea + B'^Bay/^ (16) 

C+{t+) 

We also define Q' on all C^(t~) by Lie transport of Q on C~ and define for all 

t~ < similarly. 

In the following, we shall restrict consideration to spacetimes which satisfy the following 
"decay hypothesis": 

Decay hypothesis: A spacetime (M, Qab) satisfying the conditions stated at the beginning 
of this section will be said to satisfy the decay hypothesis for Maxwell fields if for every 
smooth (in the physical spacetime) solution, Fab, of the source-free Maxwell equations which 
satisfies our Maxwell field assumptions, we have 

lim = lim ||F||(t") = (17) 

f+— »oo t~—f—00 

It should be noted that although the definitions of HFlKt"*") and depend upon 

the choice of conformal factor, Q, on and C~ , it is clear that satisfaction of the decay 
hypothesis does not depend upon this choice. Furthermore, it is not difficult to verify that 
the satisfaction of the decay hypothesis does not depend upon the choices of C"*" and C~ . 
Thus, the decay hypothesis is, indeed, a condition on the spacetime, {M,gab)- 

If the Killing field t° is spacelike somewhere in M (i.e., if an ergoregion is present), then 
solutions of the source-free Maxwell equations with negative total energy can be constructed, 
and since only positive energy can be radiated to infinity, it is clear that the decay hypothesis 
cannot hold (see |1^). However, we conjecture that the decay hypothesis holds for all 



spacetimes satisfying our spacetime assumptions in which is globally timelike. 

Note that if (M, gab) satisfies the decay hypothesis, then there can exist at most one 
stationary solution to Maxwell's equations with a stationary source ja which satisfies our 
Maxwell field assumptions. (Proof: If two such solutions existed, their difference would be 
a source-free solution which does not decay.) For a stationary ja, both the advanced and 
retarded solutions are necessarily stationary. If ja is stationary and has compact spatial sup- 
port on a Cauchy surface, then condition (1) of the Maxwell field assumptions automatically 
is satisfied. The final property we shall need is that the advanced and retarded solutions 
associated with such a ja satisfy conditions (2) and (3) as well: 

Stationary solution property: (M, gab) will be said to satisfy the stationary solution 
property if for any stationary ja of compact spatial support, the advanced and retarded 
solutions satisfy conditions (2) and (3) of the Maxwell field assumptions. 

It is readily verified that the stationary solution property holds for Minkowski spacetime, 
and we believe that the stationary solution property holds for all spacetimes satisfying our 
spacetime assumptions. Indeed, near infinity, the behavior of stationary, source free solutions 
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is such that in order to prove that this property holds, we would, in essence, need only show 
that the advanced and retarded solutions for a stationary ja of compact spatial support do 
not blow up at spatial infinity. While this seems undoubtedly true, one would need some 
bounds on the long distance behavior of the advanced and retarded Green's in a general 
spacetime satisfying our assumptions in order to obtain a proof. 

Note that, in view of the above uniqueness property, if {M,gab) satisfies the decay hy- 
pothesis and the stationary solution property, then for any stationary ja of compact spatial 
support, the advanced and retarded solutions are equal. Note also that for any stationary 
solution, the flux of energy through X+ between and any cross-section (or through I~ 
between and any cross-section S~) must either be zero or infinite. The latter possibility 
is ruled out by our Maxwell field assumptions. Thus, when the stationary solution property 
holds, the retarded (= advanced) solution for any stationary ja of compact spatial support 
must have vanishing energy flux through T"*" and X~ . From the form of the stress-energy 
tensor, it can be seen that this is equivalent to the vanishing of the pullback to X"*" (or X~) 
of FabU^, where n"- denotes the normal to X+; equivalently, we have on X~*" 

F^^n' = afia (18) 

for some function a on X+. 

We have the following theorem: 
Theorem 2.2: Let (M, gab) be a spacetime satisfying the conditions stated at the beginning 
of this section which, in addition, satisfies the decay hypothesis and stationary solution 
property. Let Fab be a smooth solution to Maxwell's equations satisfying our Maxwell field 
assumptions such that, in addition, ja— ja vanishes outside of a compact region of M, where 
jf is a stationary charge-current distribution of compact spatial support. Then, we have 

j^t^'Tabtcde' — j^t'^Tab^-cde' = J^t"' Fabj'' Gcdef , (19) 

Proof: Without loss of generality, we may assume that the support of on S is contained 
in S', so that j^ vanishes in the region, V{0, 0), bounded by and C~ . Since eq. (^ holds 
for Fab, it suffices to show that 

lim / eTabecde' 

t+^+oo Jc+{t+) 

exists and equals 

lim / t'^Tabe.de (21) 

t-^-oo Jc-{t-) 

Let FJ, be the retarded (= advanced) solution to Maxwell's equations with source j^, and 
let 

F'ab = Fab - -^afe (22) 

Then, by straightforward estimates similar to those used to obtain eq. (0) above, we obtain 
for all t+ > 0, 

9 



(20) 



IC+(t+) 



t Tab^cde 



^arpS ^ b 



C+{t+) 



- ab^cde 



< K 



\F'\\it' 



(23) 



for some constant where T^j, denotes the stress-energy tensor of Fj,. However, by sta- 
tionarity, we have 



j^arpS b 



C+{t+) 



■ ab^cde 



c+ 



ab^cde 



(24) 



Furthermore, since F' ab is source free for sufficiently large t^, the decay hypothesis applies 
to it. Consequently, we obtain 



Similarly, we obtain 



lim / t'^Tabe.de 
i+^+oo Jc+{t+) 



lim / f'Tabecde' 



c+ 



j.arpS b 



t-^-oo Jc-{t-) 



c- 



±arpS b 



(25) 



(26) 



Finally, we apply eq. (0) to Fj,, choosing t+ = t~ = 0. Taking into account the fact that 
the energy flux integrals through X"*" and X~ vanish for and that vanishes in 1^(0, 0), 
we find that the right sides of eqs. ( 



I) and (|26D are equal, as we desired to show. 
Applying this theorem to the case of a source-free Maxwell field Fah-, we have 



□ 



/ t^'Tab^cde ~ / t"'Tabecde — 0. 
Jl+ Jl- 



(27) 



That is, for a source-free solution, the energy radiated "into" the spacetime through past 
null infinity is equal to the energy radiated "out of" the spacetime through future null 
infinity. This observation, which will be important in our subsequent analysis, demonstrates 
the intuitive content of our decay hypothesis: that no energy remains "trapped" in the 
spacetime. 

We are now prepared to generalize the above theorem to the case of a point particle source 
with charge e and world line z{t). In this case, the current ja is given by the distribution 



ja{x) = e 5{x,z{T))UadT. 



(28) 



Since Maxwell's equations are linear, they are well defined for distributional sources. We 
shall consider only distributional solutions Fab which are smooth away from the world line 
of the particle and which satisfy the Maxwell field assumptions stated near the beginning 
of this section. In that case, the integrals at X"*" and X~ appearing in eq. (p!9[) are well 
defined. However, since Fab is necessarily distributional in a neighborhood of the world line 



of the particle the integral on the right hand side of eq. (19), which represents the work 
done by the source on the field, contains a formal product of distributions and is therefore 
ill defined. Nevertheless, we now wish to show that, in the point particle case, we obtain a 
new version of eq. (|19]) in which the problematic product of distributions on the right hand 
side is replaced by the integral 
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dr, (29) 



where0 



2 1 

fa = eF'^v!" + -e^{<Xa - a^Ua) + -e^{RabU^ + UaRbcU%^) 



+e^v' r 2V[aGfyU''{T')dT', (30) 

is the expression given by DeWitt and Brehme (as corrected by Hobbs) for the total electro- 
magnetic force on a point particle coupled to a Maxwell field Fab- In this force expression, 
is the retarded Green's function for the vector potential, satisfying 

V'VbG^a' - Ra'Gfa> = -^7rgaa'S{x, z) , (31) 

F^ is the retarded solution with source (|28|), and 

Kb = Pah - F^h- (32) 

In the integral over the particle's past world line, usually referred to as the "tail term" , there 
is an implicit limiting procedure: the integral is performed from — oo to a point f < r and 
then the limit f — >■ r is taken. One consequence of the Hadamard expansion is that this 
limit exists; i.e., no singular part of Gaa' is encountered in this limit. 

In order to better understand the behavior of eq. (|TU]) in the point particle limit, we 
write Fab as the sum of the half-advanced plus half retarded solution 

Fab^{Fl + Ftb)/2 (33) 

and a source-free solution 

Fir = Fab - Fab. (34) 

We shall consider only particle motions such that F^ and F^ are smooth away from the 
world line of the particle and such that F^ and F^ satisfy our Maxwell field assumptions. 
If the stationary solution property holds, this will hold whenever there exist t~, such that 
the particle is stationary for all t < t~ and all t > . 
Since Tab is quadratic in Fab, we have 

Tab[F, F] = Tab[F, F] + Tab\F'''% F^"--] + 2Tab[F''''% F], (35) 

where we have defined 



^The sign of the "tail term" in eq. (19) of reference |^ is incorrect due to an error in transcribing 
DeWitt and Brehme's expression into our notation. Additionally, a factor of 2 error was introduced 
in copying the "tail term" to eq. (23). As a result, the "tail term" in eqs. (25) and (26) of 
reference |5| must be multiplied by a factor of —2 in order to produce the correct formula (given 
in eq. (|30|) of the present paper). 
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T^,[F\ F^] = - i^Fl^%^ - -ga,{Fy'F'^^ . (36) 

Now, F^^"^ is a source-free solution which is smooth away from the world line of the par- 
ticle and consequently (by the propagation of singularities theorem [^) is smooth every- 



where. Since F^^^"^ also satisfies our Maxwell field assumptions, eq. (|27|) applies to it, and 
Tab[F^'^'^'^ , F^^*^^] will make zero contribution to the net energy radiated to infinity (the left 
hand side of eq. (|19D above). Furthermore, since TahlF^^^"^, F] is the product of a distribution 
and a smooth tensor field and thus is well defined as a distribution, this term will give 
rise to a well-defined volume integral when we apply Stokes' theorem; we will evaluate this 
term in Proposition 2.1 below. Thus, the only term that is mathematically problematic is 
Tab[F, F], which contains a product of distributions. However, we will now show that, if the 
point particle source is initially and finally stationary (i.e., if there exist such that the 

particle is stationary for all t < t~ and all t > t^), this term makes vanishing contribution 
to the net energy radiated to infinity. 
To show this, we write Tab[F, F] as 

Tab[F, F] = Tab[r^\ r^] + r,4F^, f% (37) 

where 

= {Fl - Ftb)/2. (38) 

Like F^^'^, the radiative solution Fl^ is a source-free solution which is smooth away from 
the world line of the particle and, consequently, is smooth everywhere. Since Fl^ also 
satisfies our Maxwell field assumptions, Tah[F'^^'^ , F^^] makes no contribution to the net 
energy radiated to infinity. Therefore, we have 

/ t"'Tab[F , F]ecde — I t°'Tab[F,F]ecde = / t'^Tab[F^,F^]ecdJ' — I f^TablF^ , F^Ycde 
Jl+ Jl~ Jl+ Jl- 

(39) 

However, since the particle is initially stationary, it follows that F^ satisfies eq. (0) on I~. 
Taking into account the fact that on t°- is proportional to the normal, n"-, to and F^ is 
antisymmetric in its indices, we see that the flux integral over X~ of the "advanced-retarded 
cross-term" vanishes. Similarly, using the fact that the particle is finally stationary, we see 
that the integral of this cross term over X+ also vanishes, which yields the desired result. 

We are now ready to state and prove the following intermediate result^: 
Proposition 2.1: Let {M,gab) be a spacetime satisfying the conditions stated at the be- 
ginning of this section together with the decay hypothesis and stationary solution property. 



^It should be noted that a direct analog of this proposition also holds in the case of smooth sources 
and, indeed, the present proposition could be obtained by taking the point particle limit of the 
analogous result for smooth sources. We chose not to present this result in the smooth case because 
Theorem 2.2 (i.e., the smooth source analog of Theorem 2.3 below) could be proven more directly 
by other means. 
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Let z{t) be a timelike curve which differs from an orbit, zq{t), of the stationary Kilhng field 
t°- only over a finite interval. Let Fab be a solution to Maxwell's equations with source (^Sj) 
which satisfies our Maxwell field assumptions. Then, we have 



f'Tatecde" - / t^Tate^J = / f^fadr, (40) 
x+ Ji- J 

where 

fa ^ eF^ru' (41) 

Proof: We apply Stokes' theorem to the differential of energy current three-form associated 
with the cross-term 2Tab[F^'^^'^ , F]. In analogy with the derivation of eq. (j|), we obtain 



hm 2/ eTab[F''^^,F]e,J- lim 2/ i^]ecde' 

^^+00 JC+{t+) t-^-00 Jc-{t') 

+ 2/ tX6[^'"',^]eV-2 / rTab[F''''^.F]e,J = -2 [ t'^V'Tat[F'''^^ F] 
Ji+ J I- Jm 

- t I'ab J ^cdef 
JM 

eF^^^^u^dr 
However, by the above results, we have 

Thus, it suffices to show that 

lim / ^'^^„4F^-^F]eV= lim / eTab[F'''\ F]e,J = Q (44) 

t+^+ooJc+{t+) t- -*-oo J C-{t-) 

To show this, in analogy with eq. ( ^21) we define 

F"ab = Fa, - F!, (45) 

where F^j, is the stationary solution sourced by a point charge following the orbit, Zo{t), 
of t". We use eq. (|45|) to substitute for F in eq. (^4[) . Since F^'''''^ is a source free solution 
satisfying our Maxwell field assumptions and, in the asymptotic past and future, so is F", 
it follows immediately from the decay hypothesis that 

lim / eTab[F''"',F"]e,,e'= lim / rT,4Ff-^ F"]eV = (46) 

t+^+oo Jc+it+) t-^-ooJc~{t~) 

On the other hand, although F^ is singular on the world line of the source, the Hadamard 
expansion shows that on C~^{t~^) it diverges as 1/a as one approaches the point charge, 
where a denotes squared geodesic distance. Consequently, the unphysical orthonormal frame 
components of Fj, are functions on C^(t^) with respect to the unphysical volume element 
^abcdn""- Therefore, we obtain 
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C+{t+) 



IC+{t+) 



rifree zpSl- b 
-r Cede 



Ui||^^''i (47) 
for some constant k, which, together with the decay hypothesis apphed to F^'^^'^, imphes 



lim / eTab[F 

t+^+oo Jc+(t+) 



free tjiSi ^ b 
-T Cede 







(48) 



and similarly for the limit as t ^ — oo. g 

The above proposition is a somewhat curious result, because although it provides an 
expression of the general form we are seeking, the force ( ^l]) is not the force expression given 
by DeWitt and Brehme (Eq. (|30| ) above). Indeed, we have 

2 o / . o \ 1 

— ( 

3 

1 



fa 



2 

a Ua, 



2 



/T 1 /' + 00 

2V[,G5,,n^ (r') rfr' - -e%' / 2V[,G^,,«^ (r') dr', 
-oo ^ -/r 



(49) 



which differs from eq. (|30|) by the form of the tail term. Indeed, for the retarded solution 
(so that vanishes), eq. (^) contains acausal contributions from the future history of 
the particle, whereas eq. (|30|) does not. Thus, although fa is the correct force expression 
in Minkowski spacetime (where the tail term vanishes) and it has been used in curved 
spacetime, most notably by Gal'tsov , apparently because of the calculational simplicity it 
affords, it is not equivalent to fa- Indeed, we always have fa = whenever Fal^^ = 0, so that, 
in particular, /a = for a static charge in Schwarzschild if there is no incoming radiation. 
On the other hand, it is not difficult to verify that the force prescription given by Smith and 
Will [0 for a static charge in Schwarzschild satisfies the axioms of and thus agrees with 
the DeWitt-Brehme force. The force calculated by Smith and Will is nonvanishing. 

We now calculate the difference between the work done by fa and fa- The difference 
between fa and fa is given by 
1 



fa-fa = ^e\' / 2V[aGLM'^ (t') dr' + -e\' / 2V[,Gi,^i^ (r') dr' 

Z J— CO Z Jt 



-oo 
hoo 



where 



Gab' — {Gab' + ^afe')/2 



(50) 



(51) 



and where it is understood that the singular point r' = r is omitted from the region of 
integration, i.e., one excludes an interval of size e centered at r and then takes the limit as 
e ^ 0. Therefore, the difference between the total work done by and the total work done 
by fa during the time interval [r~,r+] isR 



^Note that we proceed by considering the work done in the finite time interval [t^,t+] and only 
at the end of the calculation do we take the limit — > ±oo. We do so in order that all of the 
integrals arising in our calculations will converge suitably well to justify the interchanges of orders 
of integration. 
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/ - fa) dT = e^ / 2VlaGb]b'u' dr' 



dr 



—e 



£tGwu\^' dr' dr. (52) 



We now analyze the contributions of the two terms appearing on the right side of eq. 

Consider the first term. On account of the finite range of the r integration, we may 
interchange the orders of integration to obtain 



/ / u'Vb{t''Gab')u' dr'dr = e^ / [eGab'^T-^ ,r') - eGav{r- ,r')]v' dr' 

Jt~ J—oo J —oo 



rAf'(z(r+))-tMr(^(r-)) 



tail / 



(53) 



where A^J^^^ is the "tail part" of the time-symmetric vector potential for our point particle 
source. (Recall that the point r' = r^ is omitted from the integration in eq. ( |5BD so that 
A^f^ is finite.) Now, since the particle is being held stationary for early and late times, the 
solution must become stationary at early and late times, soQ 



lim rAf'\z{r^)) = -y*^'i(x±, 



(54) 



where y*'^'^(xi, X2) is the tail part of —f^Aa as measured at spatial position xi for a 
point particle source held stationary at spatial position X2 for all time and where x^ = 
limT-±^-i-oo x(r^). Therefore, we have for the first term 



rT~^ roo _ 

lim eM / u^Vb{t''Gab')u^ dr' dr = -eV'^^{x+,x+) +eV'^^{x-,x- 



(55) 



Since we further assume that the particle begins and ends in the same spatial position, we 
have x+ = x~ and the right side vanishes. 

On the other hand, the second term in eq. ([5^ ) can be analyzed as follows. The Green's 
function G satisfies the property that it is symmetric in its arguments,^ 

Gaa'{x,x') = Ga'aix',x). (56) 

Furthermore, the invariance of G under the time translation isometrics implies that 



^Strictly speaking, the decay hypothesis directly tells us only about the decay of the Maxwell 
field tensor, Fab- In order to obtain eq. (|5^), we must extend the decay hypothesis to include the 
assumption that in the Lorentz gauge, the vector potential of any smooth, source free solution goes 
to zero at late times (say, pointwise along every Killing orbit). 

^This fact follows from the antisymmetry of — together with the support properties of 
G^ and G^. The antisymmetry of G^ — G^ follows directly from the generalization to curved 
spacetime of the Maxwell version of lemma 3.2.1 of |15| (see also Q). 
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£tGaa' + £t'Gaa' — 0. 

Consequently, if we integrate over the same limits in r and r', we have 

J J £tGbi,'{T,T')u^u^ dr' dr = J J £t'Gb'b{T' ,t)u^u^ dr' dr 

£t'Gbb'{T,T')u^u^ dr' dr 
£tGbb' (t, t')u^u^ dr' dr 



(57) 



(58) 



Here the first equality was obtained by interchanging the dummy variables r and r', the 
second equality was obtained using eq. (0), and the last equality was obtained using eq. (^7|). 
This shows that any integral of the form (^) vanishes, provided only that the domain of 
integration is symmetric in r and r' and that the integral is suitably convergent in this 
domain to justify the interchange of orders of integration. 

Applying eq. (^) to the domain [r_,r+] x [r_,r+], we find that the second term in 
eq. (|52D becomes 



£tGbb'U^u^ dr' dr = 



£tGbb'U^u^ dr' dr + I I £tGhh'U^u^ dr' dr 



(59) 



Focusing attention on the first term, we find 



lim / £tGbb'U^u^ dr' dr = — lim 

r+— >+oo Jt^ Jt+ r+— >+oo jt^ J t' 



£t'Gbb'U u dr' dr 



lim 



lim e 



£Au' Gbb')dT' 



vJ'dT 



Xv"Gbb'{z{T),z{T+))v'dT 



lim -xu'''A^,{t^ 

r+— >+oo Z 



(60) 



where we have used the fact that t"" = xu°' for sufficiently late times. Therefore we have 



lim e 

r+— >+oo 



£tGbb'U^u^ dr dr = — - 

oo Z 



(61) 



Again, this term vanishes when x"^ = x~. 

The above result together with the previous proposition yield the main theorem of this 
section: 

Theorem 2.3 Let {M,gab) be a spacetime satisfying the conditions stated at the beginning 
of this section together with the decay hypothesis (and its extension indicated in footnote ^ 
and stationary solution property. Let z{t) be a timelike curve which differs from an orbit, 
zo{t), of the stationary Killing field only over a finite interval. Let Fab be a solution to 
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Maxwell's equations with source p8| ) which satisfies our Maxwell field assumptions. Then, 
we have 

IfT^becJ - j^ f'T^becde = J f'fa dr. (62) 



where /° is the DeWitt-Brehme force, eq. (|30[) . 

Thus, both the correct force prescription fa and the incorrect force prescription fa exhibit 
the property of global energy conservation This underscores the fact that global energy 
conservation is insufficient to determine a local expression for the force on a point particle. 

In Proposition 2.1 and Theorem 2.3, it was required that the particle begin and end 
in the same stationary state. It is interesting to consider the case where the particle is 
stationary for t < t" and again becomes stationary for t > t^, but its final position differs 
from its initial position. In this case the proof of Proposition 2.1 holds without essential 
change, and eq. (^Oj) still applies. From eqs. (p2D, (pSf), and (0), we obtain 

Ij^Tatecde' - l^t'^Tabecde' + V*"'(x+, x+) - ^V'^^^' , x") = J fa dr. (63) 

Comparison with eq. (H) strongly suggests that we identify |K*'^'^(x, x) as the renormalized 
electromagnetic self-energy, -Egeif? of a stationary point charge at position x 

E,eif(x) = |r^"(x,x) (64) 



From the analysis of Smith and Will [T^, it can be seen that for a static point charge in 
Schwarzschild spacetime, we have 

Eseif = e^M/2r^ (65) 

where M is the mass of the Schwarzschild spacetime and r is the Schwarzschild radial 
coordinate. For sufficiently slow motion (relative to the static Killing field), the static force 
associated with eq. (^) should dominate over any "damping force" associated with radiated 
energy^. Therefore, in the slow motion limit, we would expect to obtain the dominant self- 
force correction to the motion of a freely falling test charge in Schwarzschild spacetime by 
replacing the usual expression for energy, 

E = -mgabf'u^ = m(^l- i, (66) 

by 

E' = E + E,eif . (67) 



®The static force associated with eq. (p5|) is of order e^M/r^ , whereas in the slow motion, weak 
field limit, the damping force on a circular geodesic orbit is of order e^uj'^r = e^M^/^/r^/^ (see @), 
which is smaller by a factor of (M/r)^/^. 
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since E' (rather than E) should be a constant of the motion. Solving for t and plugging into 
the expression gabU°'u^ = — 1 with L = mr'^(f), we obtain 



1 2M\ / ^ m\ ^ E^^uE' Ej,, _ E'^ ^ 

- ' ' 2mr2 2 1 m 2m 2m' 



which is the equation of motion for a particle of mass m and energy E''^/2m in a one- 
dimensional potential, mlles- To lowest nontrivial order in e, Ues is given by 

1 M L2 ML^ e^ME 
Ues = - + TT-^ ^ + TT^- (69) 



2 r 2m?r'^ m?r^ 2m?r'^ 

The last term in eq. (^) represents the "self-force" correction to the motion of a charged 
particle in Schwarzschild spacetime. Interestingly, for the case of a nearly extreme Reissner- 
Nordstrom black hole, a self-energy correction of the form (^) with E^cm given by eq. (|65|) 
is of just the right nature to eliminate the possible counterexamples to cosmic censorship 
recently proposed by Hubeny |jl9 . 



How are our arguments and results modified if the particle motion is not stationary in 
the past and/or future? In our arguments above, we used stationarity in the past and future 
only to obtain the following four results: (1) To conclude (in conjunction with our station- 
ary solution property) that the advanced and retarded solutions satisfy our Maxwell field 
assumptions. (2) To conclude that the "advanced-retarded cross-term" makes no contribu- 
tion to the energy flux through X"*" and I~ (see the right side of eq. (|39|) above). (3) To 
prove eq. (|4^) . (4) To evaluate the right side of eq. (|^). Thus, eq. (|6^) will continue to 
hold in all other circumstances where the above results are valid and where the right side of 



eq. (|52D vanishes. In particular, if we consider a particle orbit which comes in from infinity 
in the asymptotic past in (nearly) geodesic motion and emerges to infinity in the future also 
in (nearly) geodesic motion, then, by calculations similar to those done in the stationary 
case, the right side of eq. (|5^ ) can be shown to be given by the limit of "tail terms" in the 
asymptotic past and future, which should vanish. However, it is somewhat more delicate as 
to whether the other results hold. For example, results (l)-(3) should hold for a point charge 
in Minkowski spacetime which undergoes exactly inertial motion in the asymptotic past and 
future. However, for the case of two point charges in Minkowski spacetime which move near 
infinity under the influence of each other's Coulomb field in the asymptotic past and future, 
property (1) does not hold since fails to smoothly extend to X~ and fails to smoothly 
extend to X+ pO| , |2T[] . Nevertheless, the "pullback" of F^ does continuously extend to X~; 



i.e., if we contract F^ into smooth two vector fields which, at X~, are tangential to X^, then 
the resulting scalar field continuously extends to X~. Since only these components of F^ 
are relevant for our energy flux calculations, it appears that, even though result (1) fails, 
eq. (B^) will hold nevertheless in this case. More generally, we expect eq. (|52|) to hold for 



particle orbits in curved spacetime which come in from infinity in the asymptotic past in 
(nearly) geodesic motion and emerge to infinity in the future also in (nearly) geodesic mo- 
tion, provided, of course, that the spacetime satisfies our spacetime assumptions, the decay 
hypothesis, and the stationary solution property. 

If the spacetime (M, gab) is axisymmetric as well as stationary, then essentially all of the 
results leading to Theorem 2.3 carry over straightforwardly if we replace the timelike Killing 
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field by the axial Killing field (j)"'. The only exception is the analysis of the "advanced- 
retarded cross-term" of eq. (|39|) above, which now no longer automatically vanishes as a 
consequence of eq. (|1^). Nevertheless, under the additional assumptions that point charge 
stationary solutions are axisymmetric on T"*" and X~ and that the vector potential of a 
source-free solution can be chosen to vanish as one approaches timelike infinity along 
and X~, this cross-term can be seen to vanish by the following argument. 

For definiteness, let us consider the integral over past null infinity. In light of our as- 
sumption that the point particle source ja differs from the stationary source only in a 
compact region of M, on we can write and as 

F^b = f!, (70) 
F^b = F!, + (71) 

where Fj, is the retarded (= advanced) solution associated with j^. Then is source-free 
outside of a compact region of M and vanishes outside of the causal past of this region. It 
follows from eq. (p!8| ) that Tab[F^,F^] radiates no angular momentum through X~, so we 
have 

2 1^ rTabiF"", F^] e,,/ = 2 ^TablF^ F'] 

= X. r {F!,Flf + F^.F^f - IgabF^F'f^) e^J. (72) 

Since 0" lies in X~ , the last term clearly vanishes. Writing out and Fj, in terms of the 
corresponding vector potentials, we arrive at 

2 /_ rT,4F^, F^] = I (rV[,4F,'^ + ^Vi^^jF^/) e,,/ 

= Vf{rAl)F"f e,deb - £^A)F''f e^^b 

+ VfirK)iF''f^ e,deb - £c^A'^iF^ff e,aeb (73) 

The second integral is immediately seen to vanish under our assumption that the vector 
potential of a stationary solution can be chosen to be axisymmetric at . Integrating by 
parts in the fourth integral gives us 

£^A'f{F''f^ e,deb = -l_ A'f£^{Fy^ e^^eb, (74) 

which vanishes under the same assumption. (The "boundary term" makes no contribu- 
tion since the orbits of 0" are closed). For the first integral, since F'"* satisfies Maxwell's 
equations, we have 

Vfi^ADF"^ ebcde = 3 [(0 ■ A^) * F'] , (75) 

where the two-form *F' is the Hodge dual, defined by *F^^ = — (l/2)ea6'^'^F^'^. Using Stokes' 
theorem, we can convert this integral to two boundary integrals, one near and the other 
near i~ . Since the two- form *F' is source-free outside of a compact region of M, the integrand 
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vanishes near z° by causality and vanishes near i under the decay hypothesis. Similarly, 
for the third integral we have 

V/(r ■ K){Fyf ee,e6 = 3j^_d[{<p- A') * F^] . (76) 

This term vanishes because the vector potential A'^ can be chosen to vanish near by 
causality and, by our assumption, can be chosen to vanish near i~ . 

By reversing the roles of and F^ and repeating the above arguments, we also show 
that the integral over X"*" vanishes. Thus, if (M, gab) is axisymmetric as well as stationary 
and the additional conditions given above are met, then under the hypotheses of Theorem 
2.3, eq. (|6^ continues to hold if we replace by 0°. 

Finally, we note that our results can be generalized to allow for the presence of a black 
hole and white hole whose horizons, Ti^ and 7i~, intersect on a compact 2-surface, S. We 
again require (M, gat) be a globally hyperbolic, stationary spacetime that is asymptotically 
fiat at null and spatial infinity. However, we now require that there exist a spacelike hyper- 
surface, S, with boundary, S, of the form S = Send U S' with S' compact, such that S \ S* is 
a Cauchy surface for the domain of outer communications. In this case, we deform S to the 
future in small neighborhoods of both i° and 5* so that the deformed surface, C~^, intersects 
both X"*" and the future event horizon H"*" in a cross-section. We construct C~ similarly. 
With these changes, our previous analysis for the case of no black or white holes can now 
be repeated^ wherein in all of our assumptions, arguments, and formulas, we replace X+ 
by X+ U Ti^ and X^ by X^ U Ti.^ . In particular. Theorem 2.3 holds for the above class of 
spacetimes containing black holes and white holes provided only that we modify eq. ( |62l) so 
as to include contributions from the energy flux through 7i"*" and 7i~ in exact parallel with 
the contributions from X"*" and X~. 

However, it would appear that some significant differences do occur in the classes of 
spacetimes for which the decay hypothesis and stationary solution property hold. Recall 
that when no black hole or white hole is present, we conjectured that the decay hypothesis 
would hold if and only if no ergoregion is present. However, for a rotating black hole — i.e., 
if is not normal to Ti^ and 7i~ — the "ergoregion instability" argument [l^ no longer 



applies, since negative energy can be radiated into the black hole. Thus, it seems plausible 
that the decay hypothesis will hold for many rotating black holes — in particular, for Kerr 
black holes — despite the presence of ergoregions. On the other hand, for the case of a 
rotating black hole, it does not seem plausible that the stationary solution property holds 
(where, with the modifications described above, this property now includes the requirement 
that the Maxwell field extend continuously to 7i+ and H^). Specifically, for the retarded 
solution for a non-axisymmetric, stationary source, there should be a steady, nonvanishing 
flux of angular momentum into the black hole. Consequently, F^ should not be continuously 



^The only exception is that for a rotating black hole, the vanishing of the advanced-retarded 
cross-terms does not follow from the horizon analog of eq. (^) but instead must be shown by the 
arguments similar to those used above for the case of angular momentum. However, as we shall 
note shortly, the validity of the stationary solution hypothesis appears to be limited to non-rotating 
black holes in any case. 



20 



extendible to 7i~. Similarly, should fail to be continuously extendible to TC^. However, 
it seems plausible that the stationary solution property always holds for a non-rotating black 
holef^. Thus, it appears likely that the hypotheses of Theorem 2.3 hold precisely for the 
case of non-rotating black holes that contain no ergoregion. Of course, the conclusion of 
that theorem may still hold under more general hypotheses. 

Finally, in the axisymmetric case, we may consider conservation of angular momentum 
in spacetimes containing a non-rotating black hole and white hole. As in the case where no 
black hole or white hole is present, all of our arguments carry through straightforwardly when 
t"' is replaced by 0" except for the vanishing of the advanced-retarded cross term. However, 
the presence of a black and white hole now produces a significant difference because, even 
in Schwarzschild spacetime, the field of a stationary point charge fails to be axisymmetric 
on the horizon. Thus, the symmetry argument used above to obtain the vanishing of the 
advanced-retarded cross terms at T"*" and X~ is inapplicable to and H~ . We have not 
been able to show that these cross terms vanish (or that the cross term contribution from 
7^+ cancels that from Ti^). It would be of interest to determine if any violations of angular 
momentum conservation can occur in this case and, if so, how large the violations could be. 

III. GRAVITATIONAL CASE 

In this section, we shall give an analysis of conservation results in the gravitational case 
in close parallel with our analysis of the electromagnetic case given in the previous section. 
Although there are many similarities between the gravitational and electromagnetic cases, 
there also are some key differences as well as some additional technical difficulties occurring 
in the gravitational case. Consequently, with the exception of Theorem 3.1 below, we shall 
not attempt to prove theorems based upon analogs of our Maxwell field and decay hypothesis 
assumptions, but rather will merely sketch how it should be possible to obtain similar results. 
In this sense, our results for point particles in the gravitational case will be considerably 
weaker than in the electromagnetic case. 

To begin, we wish to analyze energy conservation for solutions to the linearized Einstein 
equation on a spacetime, {M,gab), which satisfies the spacetime assumptions stated at the 
beginning of the previous section. However, unlike the electromagnetic case, in order to have 
well defined linearized field equations, it is necessary here that the background spacetime be 
a solution to Einstein's equation. Although there should not be any difficulty in principle 
in allowing the background spacetime to possess suitable matter fields, in order to directly 
apply the results of in the case of point particles, we must restrict to the case where 
the spacetime is vacuum at least in a neighborhood of the worldline of the particle. This 



^'^In particular, in the static case, the equality of and F^ for static sources throughout the 
domain of outer communications follows from the time reflection symmetry. Since F^ is regular 
on Ti.'^ and F^ is regular on TC^ it follows that both must have continuous limits to 7{~^ U 'H~ , in 
accord with the stationary solution property. Note, however, that F^ will be discontinuous across 
Tl~ and F^ will be discontinuous across Ti'^ in a manner similar to the behavior of the retarded 
and advanced solutions for a uniformly accelerating charge in Minkowski spacetime [p^. 
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restriction would still allow us to consider, for example, a fluid star solution in which the 
particle remains exterior to the star, but the inclusion of fluid or other matter flelds would 
signiflcantly complicate the analysis. Hence, for simplicity, we shall restrict attention to the 
case where the background spacetime satisfles Gab = 0. Unfortunately, this restricts the 
background spacetime to be Minkowski spacetime.Q More generally, when black holes are 
admitted, the stationary black hole uniqueness theorems restrict the background spacetime 
to being Kerr. However, we shall not make explicit use of any special properties of the back- 
ground spacetime in our analysis below, and we believe that our analysis should generalize 
to cases where suitable matter is present in the background spacetime. 

On (M, Qab), we initially wish to consider smooth solutions, jab, to the linearized Einstein 
equation with a smooth source Tab 

G'l^M = 87rr„, (77) 

Our basic strategy here will be to repeat the analysis of the previous section, replacing 
the electromagnetic energy current three-form f^Tab^cde with the "effective gravitational 
energy current" three-form t^'Tab^cde' = ~^t°'G'^ab^cde' constructed from the second order 
Einstein tensor G^^-* (see [|lOl) associated with the linearized solution. In doing so, we will 
rely on the results of Habisohn , which show that the net flux from this effective energy 
current at null infinity agrees with the net Bondi flux to second order. However, in order 
to apply Habisohn's results, we need to assume that Tab vanishes in a neighborhood, U , of 
X"*" U T~ U z° that contains complete Killing orbits and that, in [/, 7^6 satisfies both the 
weak and strengthened fall-off conditions of [|IU|. Here, the weak fall-off conditions of M 



require that vanish in a neighborhood of and satisfy the decay properties in M fl ?7 
stated in The strengthened fall-off conditions require that, when transformed to the 
Geroch-Xanthopoulos gauge the unphysical metric perturbation 7ab = f2^7afe (which 



smoothly extends to X"*" and X~ in this gauge) satisfies similar decay properties in M fl t/, 
where M denotes the conformally completed spacetime. The requirement that 7af, vanish 
in a neighborhood of is much stronger than condition (3) imposed in the electromagnetic 
case, and, indeed, is too strong to enable us to use it directly for our purposes below. 
Undoubtedly, Habisohn's conditions could be weakened considerably, but this would require 
further analysis that we do not wish to undertake here. 

Let QabW = Qab + ^lab, whcre Gab[gcd\ = and 7ab satisfies eq. (^). The Bianchi identity 
to second order in A yields 

[vi^)G"^^(A)] = G(2)»n7c.] + V«G(i)'^n7c.] = 0; (78) 

i.e., 

V,G(2)'^^[7,,] = -C\dbcdW^'''bcd] - C\dbcd]G^'^'"'bcd] (79) 



This follows immediately from the positive energy theorem together with the fact that the 
(Komar) mass of a vacuum spacetime satisfying the spacetime assumptions stated at the beginning 
of Sec. II is easily seen to vanish. 
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where 



C\c bed] = \9'"' [Vblcd + Ve7fed - Vdlbc] (80) 

and where V a denotes the derivative operator associated with Qab- Since satisfies the 
hnearized Einstein equation ([77D, we find 

^^Q{2)ab ^ _^-^^C\,T'' - C^,r"=] (81) 

Equation ( pT] ) is a direct analog of eq. in the electromagnetic case. Furthermore, 
although in this case the physical "effective energy current three-form" 

t^Tab^cde' = — —t"" G'^^b ^cde (82) 



need not continuously extend to null infinity, Habisohn |]T0[ has shown that if his weak 
and strengthened fall-off conditions hold, its puUback to certain timelike 3-surfaces has a 
continuous limit as these surfaces approach null infinity. Consequently, in exact parallel 
with the derivation of eq. (P), we obtain the following theorem: 

Theorem 3.1 Let (M, gab) be a vacuum solution of Einstein's equation satisfying the space- 
time assumptions^ of Sec. II. Let be a smooth solution to the linearized Einstein equa- 
tion with smooth source Tab such that Tab vanishes in a neighborhood, f/, of T^ U I~ U z° 
that contains complete Killing orbits and such that, in [/, 7^^ satisfies both the weak and 
strengthened fall-off conditions of 0]. Then, 



lim / t^-Tabecde - lim / f'Tabecde 

t+^+ooJC+{t+) t-~*~oo Jc-{t-) 

+ j^t'^Tab^cde — j^t'^'Tab^cde = — J^ta[C^ bcT'^"' + C"" bcT^^]edefg- (83) 

provided that the integral over M and the limits of the integrals over C+(t+) and C^{t^) 
exist. Here, the integrals over X"^ and I~ are understood to mean the limits of the integrals 
over suitable timelike surfaces in the physical spacetime which approach and X~; these 
limits were proven by Habisohn |]lO| to exist. 

The second pair of terms in eq. ( ^5] ) was proven by Habisohn [|T^ to agree, to second order 
in 7af,, with the net Bondi energy radiated to infinity. In analogy with the electromagnetic 
case, it would be natural to attempt to interpret the first pair of terms as representing 
the difference between the initial and final gravitational energy of the matter distribution, 
and, similarly, to interpret the right side as representing the work done by the matter 
on the gravitational field. However, unlike the electromagnetic case, these terms are, in 
general, gauge dependent. Thus, unless some further restrictions and/or gauge conditions 



^^As previously noted, these assumptions actually imply that (M, g^fe) is Minkowski spacetime. 
However, as also previously noted, the theorem can be generalized to allow for black holes (in which 
case the Kerr family of solutions would be admissible) and, presumably, also could be generalized 
to allow solutions to the nonvacuum Einstein equation with suitable forms of matter. 



23 



are imposed, such an interpretation would not, in general, be meaningful. Nevertheless, 
under the hypotheses of the theorem, the second pair of terms is gauge independent — 
provided, of course that one stays in a gauge where the weak fall-off conditions hold. 

In order to proceed further, in analogy with the analysis given for the electromagnetic 
case, we need to restrict consideration to situations where the first pair of terms in eq. ( pSf ) 
cancel. As in the electromagnetic case, we consider the case where Tab differs from a station- 
ary matter distribution only in a compact region of spacetime. However, in order that '-fab 
satisfy the decay properties needed for cancellation of the first pair of terms in eq. (0), it 
also is necessary to impose gauge conditions on •jab- Note that if a choice of gauge compatible 



with the Habisohn fall-off conditions |T0| has been made so that the first pair of terms in 
eq. (|83|) cancel, the right side of that equation automatically must be gauge invariant with 
respect to any remaining restricted gauge freedom. Therefore, in that case it would seem 
appropriate to interpret the right side of that equation as corresponding to the work done by 
the matter on the gravitational perturbation, 7,16, from the "viewpoint" of the background 
spacetime, {M,gab). 

A choice of gauge which appears to be suitable for the above purpose is the Lorentz 
gauge condition, 

V^Tab = (84) 

where 

%b = {lab - ^igab) (85) 

In the Lorentz gauge, ■jab satisfies a wave equation, 

Wc%b - 2R\b%d = -l6nTab (86) 

and it seems highly plausible that, in a wide class of spacetimes, suitable decay properties 
at late and early times will hold for source free solutions with suitable fall-off near i^. In 
particular, it seems plausible that for source free solutions with suitable fall-off near i^, in 
this gauge the supremum over C~^(t~^) of the components of jab and its first two derivatives 
will go to zero as t+ — 00, and that similar decay properties will hold as ^ — 00. It also 
seems plausible that, for such solutions in this gauge, the stronger fall-off conditions of fl^ 
near T"*" and X~ will hold as —>■ 00, and t~ —00. Finally, it appears highly plausible 
that, in the Lorentz gauge, an analog of the stationary solution property (see Sec. II) will 
hold. Consequently, a direct analog of theorem 2.2 should hold in the gravitational case 
when we impose the Lorentz gauge condition. However, we shall not attempt to state and 
prove such a theorem here, since it would require considerable further analysis to make a 
judicious choice of mathematically precise formulations of the gravitational versions of the 
decay hypothesis and the other relevant assumptions needed for the proof of such a theorem. 

We may now attempt ask in the gravitational case whether, for a point particle source 
whose motion differs from a Killing orbit for only a finite time interval, the analog of eq. (|62D 
holds, i.e., whether 

'^^frabe'cde - j^t''rabe\de = J f'fa dr. (87) 
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where the integrals over X"*" and T have the same meaning as in theorem 3.1, and where 
now denotes the gravitational self- force 0,^, 

-^rn^u^'u'' £^ V^GL'fc' - ^h{G'^)ca'y) u''' u'' dr' (88) 

Unfortunately, however, this question as posed has only a trivial domain of applicability: 
As emphasized in [|^, in order for a solution to the linearized Einstein equation with a point 
particle source to exist, it is necessary that the particle move on a geodesic. Therefore, a 
point particle cannot follow a Killing orbit for a finite time interval unless that Killing orbit 
is a geodesic, in which case the particle's motion can never deviate from the Killing orbit. In 
other words, the constraints on the motion of a particle imposed by the linearized Einstein 
equation preclude the possibility of having a particle which is nonstationary for only a finite 
time interval. Note that this contrasts sharply with the situation in the electromagnetic 
case, where Maxwell's equations impose no constraints on the motion of a charged particle. 

Nevertheless, it is worth noting that it is not manifestly inconsistent to have general 
distributional sources containing point masses which are stationary outside of a compact 
spacetime region. In other words, we may have a system of point masses connected by 
devices such as rods, springs, or strings such that the point masses are initially and finally 
in the same stationary position, but are set into motion at intermediate times. Under decay 
assumptions of the type expected to hold in the Lorentz gauge as described above, it then 
should be possible to prove a direct analog of Proposition 2.1. Namely, if T^b now denotes a 
distributional stress-energy tensor which is stationary outside of a compact spacetime region, 
and if 7a{, denotes a (distributional) solution to the linearized Einstein equation with source 
Tab such that is smooth near infinity and satisfies suitable decay properties near then 
it should follow that 

J^t"'Tab<^cdJ' - j^t"'Tab(^cde = " j J a[{,C^^'"'f bcT''"' + {C^^'''')"' bcT^''](^def g ■ (89) 

where {C^'^'^^Ybc is given by eq. (pO|), with replaced by 7ab'*^, where 

7fr = 7a6-7a6 (90) 

with 

7a.= ^[7.'^, + 7a^.]. (91) 

Note that in eq. (|89D it is, in general, necessary to take into account the effect of the rods 
and other devices both directly with regard to their contribution to Tab and indirectly with 
regard to their contribution to the gravitational perturbation ^ab- 

Although we cannot even pose a meaningful question for a system composed purely of 
point masses which are stationary in the past and future, we can consider energy conservation 
for a freely falling point mass which follows a trajectory which starts at infinity and returns 
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to infinity. In a general, curved spacetimef^ such a particle of mass m will radiate energy 
of order to infinity, and one may ask if this agrees, to lowest order, with minus the net 
work done on the particle by the gravitational self-force. On account of this gravitational 
self-force, the particle will fail to move on a geodesic, resulting in self-consistency issues with 
regard to treating its motion that were discussed in However, in order to calculate the 
work done by the gravitational self-force to order m^, the deviations from geodesic motion 
produced by the gravitational self-force can be neglected. Thus, we may ask whether eq. ( ^7\ ) 
holds for a particle moving on a geodesic. 

In parallel with our discussion in the electromagnetic case (see the paragraph below 
eq. (p9| ) of Sec. II), if suitable decay properties at late and early times hold for source free 
solutions, if the advanced and retarded solutions for the particle moving on the given geodesic 
are suitably well behaved, and if the "advanced-retarded cross-term" makes no contribution 
to the flux of t°'Tab€cde' through T"*" and X~, then eq. ( ^9]) should hold, i.e., 

j^fTabtcde - j^t'^Tatecde = J fa dr. (92) 

where 

+ ^mW'^l^^ (^V'GL't' - V,(G^)e"a'6') W^'u'' dr' 
However, the difference between and is given by 

fa _ fa ^ ^2^b^c y'+J (Iv'Ghca'b' " Vfe^^'fe') dr' (94) 

Our task is to show that, just as in the electromagnetic case, this difference between the two 
forces does no net work over the world line of the particle. 
In analogy with eq. (|52D , we have 

rt''{fa-fa)dT 



(93) 



m 



m 



-00 

/r+ i'+oo / ][ _ _ _ \ ^ ^ 

- J V^^i^bcb'c' - Gacb'c'Vbt'' - fVbGcab'c'j u^u^u^ dr' dr 



00 / ]^ _ _ 

aGbcb'd — ^bGcab'c' I U u'^ dr 



u^u^dr 



m 

T Jr- 



+00 



£iGhrb'riU^u^u^ dr' dr 



00 
+00 



-m 



/ / u^Vbit''u^Gacb'd)u''u'' dr'dr 



(95) 



^■^As already mentioned, our above spacetime assumptions restrict us to flat spacetime, but the 
remarks here are relevant for generalizations to allow for the presence of black holes as well as 
generalizations to non- vacuum cases. 
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By arguments similar to those which produced eqs. (^) and (0), both of these integrals 
can be reduced to expressions involving "tail contributions" on the worldline of the particle 
for r > and r < r~. Since here we restrict consideration to geodesic particle trajectories 
which start at infinity and return to infinity in asymptotically flat spacetimes, these tail 
contributions should vanish in the limit ±00. Thus, eq. ( ^Tf ) should hold in this case, 

although the reader surely will have noted the numerous gaps in our arguments that would 
have to be filled in order to convert our arguments into a theorem. 

Finally, we comment that, in parallel with the electromagnetic case, our arguments 
should generalize to yield a conservation of angular momentum result in the case of a space- 
time that is axisymmetric as well as stationary. Similarly, our arguments should generalize 
straightforwardly to the case where black holes and white holes are present provided that the 
fluxes through the black hole and white hole horizons are included in the manner described 
near the end of Sec. II. 
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APPENDIX A: MASSLESS KLEIN-GORDON SCALAR CASE 

In this appendix, we outline the derivation of conservation results for the case of a 
massless Klein-Gordon field, satisfying 

V^Vac/. = -J. (Al) 

To a large extent, the analysis of energy and angular momentum conservation for this case 
closely parallels the electromagnetic case analyzed in Sec. II, although there are some com- 
plications arising from the lack of conformal invariance of the theory and the important 
difference that the Klein-Gordon force need not be orthogonal to the four-velocity. 
The stress-energy tensor of the massless Klein- Gordon scalar field is given by 

Tab = V,0Vb0 - ^gab9'''Vc(pVd(p (A2) 
By virtue of eq. ( [All) , it satisfies 

V'Tab = -3Va(j). (A3) 

In precise analogy with the electromagnetic case, we define the energy three-current three- 
form by 

Jcde = i°'Tabecde'- (A4) 

It follows that 

{dJ)cdef = 3t''Va(l)ecdef- (A5) 
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The task of determining falloff conditions on at null infinity is less trivial here as 



compared with the Maxwell case since the field equation (^) is not conformally invariant. 
However, the standard conformally invariant modification of eq. (|Al|), namely V^Va^ — 



= —j, is equivalent to eq. (|A1|) in vacuum spacetimes. This suggests that we define 



an "unphysical Klein- Gordon scalar field" 
conformally invariant case, i.e.. 



by using the same conformal weight as in the 



(A6) 



and require that (f) extend smoothly to X+ and I~ . We adopt this condition as the analog 
of the Maxwell field assumption (2) of Sec. II. 

In terms of unphysical variables, the physical energy three-current takes the form 



J, 



cde 



-gabg'^V cct^y d(l> 



g^^ecdef 



g'ficdef (A7) 



where Va denotes the derivative operator associated with the unphysical metric cjab (although 
it should be noted that in this formula, only derivatives of scalar fields are taken, so there 
actually is no distinction between Va and Va). Since Ua = 7^ on X = X+ UX~ we see 
that Jcde fails to extend continuously to X even when extends smoothly to X. However, 
the pull-back, Jcde, of Jcde to a surface of constant fl is given by 



J. 



cde 



^ ^ ^cde, 



(A8) 



where the functions f = ^g"'''nani, and h = VL ^f^Ua extend smoothly to X, and '^^^eabc is 
defined (up to a three- form with vanishing puUback) by 



^abcd 



n 



(A9) 



Consequently, although J^de fails to extend continuously to X, it can be seen by inspection of 
eq. (|A8|) that Jcde extends smoothly to X. Equivalently, if we choose any three smooth vector 
fields in the unphysical spacetime which are everywhere tangent to surfaces of constant Vt 
and we contract these vector fields with the free indices of Jcde-, then the resulting function 
will extend smoothly to X. 

However, it is not difficult to verify that, in general, the limit of Jcde to X will depend 
upon the choice of conformal factor, f2, used to define the surfaces to which Jcde is pulled 
back. Nevertheless, we argue now that the integral over X+ or X~ of Jcde cannot depend 
upon the choice of f2, provided only that the scalar source, j, vanishes in a neighborhood of 
X U z° and that has suitable fall-off properties at and at timelike infinity. Namely, let 
Vt and Vt' be two choices of conformal factor. Integrate {dJ)cdef over the region bounded by 
a surface of constant Q and a surface of constant Q', with suitable "endcaps" inserted near 



i and future timelike infinity. Now apply Stokes' theorem — using eq. 



the vanishing 



of j near X, and the decay of near i and timelike infinity — to conclude that, in the limit 
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as one approaches T"*", the integral of Jde over a surface of constant equals the integral of 
Jcde over a surface of constant From this it follows that the integral of Jcde over X+ is 
independent of the choice of VL. Similar results apply, of course, to integrals over T~ . Thus, 
if is smooth at X and satisfies suitable fall-off properties at i° and at timelike infinity, and 
if j vanishes in a neighborhood of X U then the total energy flux through I'^ and X~ is 
well defined. 

For a spacetime satisfying the assumptions listed at the beginning of Sec. II, we integrate 
eq. (|A5| ) over the compact region with "endcaps" given by C'^it^) and C~{t~) but now 
bounded on the "sides" by a surface of constant, nonzero VL (rather than by X U z"). We 
apply Stokes' theorem, then take the limit as ±00, and, finally take the limit as ^ 0. 
The result is the analog of Theorem 2.1, namely 

lim / rTabe^,'- lim / t^^ecde' 
Jc+{t+) t-^-00 Jc-{t-) 

+ / f'Tabecde - I t^Tabecde' = f jf'V a<pecdef . (AlO) 
JI+ JJ- JM 

where the integrals over X"*", X~, C^{t^), and C~{t~) are defined by the limiting procedure 
described above. Assuming suitable analogs of the decay hypothesis and stationary solution 
property of Sec. II, we may then obtain an analog of Theorem 2.2. 

The analog of the Lorentz force on a particle of scalar charge q in an external Klein- 
Gordon field is 

fa = q^a^ (All) 

(This equation can be derived by integrating eq. ( |A3|) over a small body and neglecting 
the self-field of the body; the Lorentz force in electromagnet ism can be similarly derived by 
integration of eq. (^).) It should be noted that, in sharp constrast with the electromagnetic 
case, the scalar force (|A11|) fails, in general, to be perpendicular to the 4- velocity, m", of the 



particle. Consequently, the rest mass, m, of the particle will vary with time, i.e., in general, 
the particle will necessarily gain or lose rest mass as a result of its interactions with the 
scalar field. Indeed, we have 

-jp = -nVa = -gM"V,0 (A12) 

Thus, the 4-momentum of a particle of scalar charge q in an external Klein-Gordon field (f) 
is given hyf^ 

Pa = {mo - q(j))ua (A13) 

where mo is a constant. 



It also should be noted that — in sharp contrast with the electromagnetic case — conservation 
of charge is not required for consistency of the Klein-Gordon equation, so there is no obstacle to 
allowing q to vary with time as well. However, we shall assume throughout this Appendix that q 
is constant. 
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When self- field effects are taken into account, the formula for the total force on a test 
particle (analogous to the DeWitt Brehme formula in the electromagnetic case) is 



fa = q^a(, 



( ^da - -(^^Ua + -RabU" + -UaRbcU^u" - —RUa 



.3 3 



6 



6 



12 



(A14) 



We wish to establish the analog of Theorem 2.3 for the scalar case, i.e., we wish to show 
that if z{t) is a timelike curve which differs from an orbit, zo{t), of the stationary Killing 
field only over a finite interval, then 



j^t"'Tab^cde' — t"'Tab€-cde — j t"" fadr. 



(A15) 



As in the electromagnetic case, we proceed by first establishing the analog of Proposition 
2.1, namely that 



JX+ 



where 



fa = q^a* 



-cde 



I free 



f'fadr 



(A16) 



(A17) 



The only significant difference from the electromagnetic case occurring here is in the analysis 
of the vanishing of the "advanced-retarded cross-term" in the integrals representing the 
energy radiated through I. This cross-term vanished trivially in the electromagnetic case, 
but here is given by a nontrivial expression. However, we can simplify this expression by 
choosing Q so that £t^ = f^Ua = 0. For the integral of the cross term over X"*", we have 



J. 



cde I 



+ /0V Va0^ + /0VVa0^ 



^ ^ (-cde 



^ ^cde 



1+ 



+ 



lA TR 



= 0. (A18) 

Here, we have used the fact that (p^ is a stationary solution near so that £t(p^ = 
f^V a4'^ = 0. In the last line, we have assumed that suitable falloff conditions on 0^ near i° 
and timelike infinity lead to the vanishing of the boundary terms. Therefore, we see that the 
"advanced-retarded cross-term" makes no contribution to the energy radiated through X. 
The remainder of the derivation of eq. (|A16|) follows in close parallel with the electromagnetic 
case. 

To compare the work done by fa with that done by we note that fa is given by 



fa = q'^a4>''' + q^ ( 7;da " ^O^a + ^RabU^ + ^UaRbcU^u'' - ^RUa 

\3 3 D D 12 



J — oo 2, Jt 



(A19) 
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Hence, the difference between the total work done by fa and the total work done by fa 
between r = and r = r"*" is given by 



t^{fa-fa)dT = q' / e 



VaGdr' 



di 



gM / eVaGdr'dr + q^ / aG dr' dr 

-q^ / eVa'Gdr'dr-q^ / fVa'Gdr'dr (A20) 



where the scalar analogs of eqs. (|57D and (|58D were used. If we now take the limit — > ±00, 
then = x(x"'")m" in the r' region of integration of the first integral, where x+ denotes 
the final spatial position of the particle and x = {~'t°"taY^'^ ■ Thus, the integrand is a total 
derivative in r'. Assuming that G falls off so that G{z{t), z{±oo)) = 0, we therefore have 



lim / e{fa - fa) dr = \im [x{^-^)G{z{t), z{r+)) - x{^-)G{z{t), z{t-))] dr 

= hm ^ rlx(x+)G^^(r),.(r+))-x(x-)G^(^(r),.(r-))]rfr, 



[x(x+)0*-i(x+,x+) - x(x-)0*^'^(x-,x-)] 



(A21) 



where ^^^^'^xi, X2) is the tail part of (p as measured at spatial position xi for a point particle 
source held stationary at spatial position X2 for all time. Thus, if the particle is initially and 
finally stationary at the same spatial position, then the work done by fa is equal to the work 
done by /„, which is in turn equal to the energy radiated to infinity by the Klein-Gordon 
field. 

As in the electromagnetic case, eq. (|A21|) suggests that for a particle held stationary at 
position X, the renormalized self-energy stored in the scalar field is given by 



E,elf(x) = |x(x)0*^'Xx,x) 



(A22) 



However, in the scalar case, additional energy also is stored in the mass of the particle. 
Indeed, from eq. ([A14|) , we obtain 



dr ^' 



a 
in 



-q- 



-q 



dr 12 



R-q^ 



idr 



G^fr, t') dr' 



d 



dr 12^ " rf^ loo^''^^' ^'^ + 



(A23) 



where in the last line, G^{t, t) should be understood to mean the limit as e — of 
G^{t,t — e). The Hadamard analysis of shows that G^(r, r) = R/12, so we obtain 



dm d . ,■ . ■,, 



(A24) 
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Thus, the mass of the particle is now given by 



m = mo - g(0^'' + 0*^") (A25) 
and, hence, the total energy of a stationary particle is 

E' = -mt^Ua + ^self 

= [mo - + + itself 

= E-|x(x)0*-i(x,x) (A26) 

where E = (mo — q<p^'^)x is the energy that the particle would have had in the absence of 
self-force effects. As in the derivation of eq. (|69|) in the electromagnetic case, in general this 
would lead to a corresponding self-energy correction to the effective potential. However, 
in Schwarzschild spacetime, Wiseman p4[ has shown that 0*^''(x, x) actually vanishes for 
all X, so, to order g^, there is no scalar self-energy correction to the motion of particle in 
Schwarzschild spacetime. 

Finally, we note that the analysis of this Appendix can be extended straightforwardly 
(in the manner explained in Sec. II) to treat the case of angular momentum and to allow 
for the presence of black holes. 
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FIGURES 

FIG. 1. The unphysical spacetime {M,gab) with the hypersurfaces C"^(t'*") and C~{t~) used 
in our analysis. 
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